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I.  Synopsis 


Research  has  uuniiiiui: d  toward  the  development  of  electro-hydrodynamically  (EHD) 
driven  micro  heat  pipe  arrays.  This  research  has  focused  on  the  following: 

•  Completion  of  the  experimentation  and  analysis  of  EHD  augmented  micro  heat  pipe 
arrays  using  dielectrophoretic  forces. 

•  Development  and  analysis  of  ion-drag  pumped  micro  heat  pipe  arrays. 

•  Demonstrated  active  thermal  control  with  ion-drag  pumped  micro  heat  pipe  arrays  for 
cyclic  heat  loadings  in  excess  of  100  W/cm2. 

The  research  has  demonstrated  that  the  forces  induced  via  the  ion-drag  mechanism  can  be 
appreciably  greater  than  by  the  dielectrophoretic  effect,  i.e.,  force  induced  as  a  result  of 
an  electric  field  present  at  a  liquid-vapor  interface.  Moreover,  because  the  ion-drag  force 
can  be  generated  within  the  heat  pipe  fluid  at  any  location,  the  localized  electric  field  can 
be  better  isolated  from  the  heat  source.  In  the  case  of  electronics  cooling,  this  is 
particularly  critical. 


II.  Primary  Goal 

The  primary  objective  of  the  research  has  been  to  achieve  active  thermal  control  of  a 
heated  substrate  for  heat  loadings  in  excess  of  10  W/cm2  using  electrostatically  induced 
forces  within  a  micro  heat  pipe  array. 


III.  Objectives 

To  realize  this  goal,  our  research  had  the  following  objectives: 

•  Investigate  the  possibility  of  using  electric  fields  to  augment  the  heat  transport 
capability  of  a  micro  heat  pipe  device  relying  upon  the  dielectrophoretic  effect  — 
e.g.,  the  force  induced  by  applying  an  electric  field  across  a  liquid-vapor  interface.  An 
analytical  prediction  of  the  transport  capability  of  this  type  of  electro- 
hydrodynamically  (EHD)  pumped  micro  heat  pipe  device  was  to  be  completed. 

•  Investigate  the  possibility  of  using  electric  fields  to  augment  the  heat  transport 
capability  of  a  micro  heat  pipe  device  relying  upon  ion-drag  pumping.  An  analytical 
prediction  of  the  transport  capability  of  this  type  of  EHD  pumped  micro  heat  pipe 
device  was  to  be  completed. 

•  Choose  the  best  of  these  two  alternative  designs  and  determine  its  potential  for 
actively  controlling  the  temperatures  of  substrates  subjected  to  cyclical  heat  loadings. 


IV.  Work  Performed 

The  following  sections  describe  the  research  that  has  been  completed  toward  these 
objectives.  The  first  details  the  accomplishments  made  in  the  development  of 
electrohydrodynamically  augmented  micro  heat  pipe  arrays.  The  second  details  the 
accomplishments  made  in  the  development  of  ion  drag  assisted  micro  heat  pipes  and  the 
success  achieved  toward  using  these  for  active  thermal  control  for  variable  heat  load. 


Experimental  and  Analytical  Studies  of  the  Maximum  Heat 
Transport  Capacity  in  Electrohydrodynamically  Enhanced  Micro 
Heat  Pipes 

Abstract 

An  experimental  investigation  was  conducted  to  evaluate  the  potential  benefits  of 
electrohydrodynamic  (EHD)  forces  on  the  operation  of  micro  heat  pipes.  In  these 
experiments,  electric  fields  were  used  to  orient  and  guide  the  flow  of  the  dielectric  liquid 
within  the  micro  heat  pipe  from  condenser  to  evaporator.  The  experiments  indicate  the 
heat  transport  capability  of  the  EHD  micro  heat  pipes  is  increased  by  up  to  6  times  of  that 
of  conventional  ones.  In  parallel,  an  analytical  model  was  developed  to  predict  the 
maximum  heat  transport  capability  for  various  electric  field  intensities  and  micro  heat 
pipe  geometries.  The  analytical  model  agrees  well  with  the  experimental  results  for  the 
geometry  studied  experimentally.  The  model  shows  that  large  pore  sizes  are  optimal  from 
a  heat  transport  capacity  perspective.  Finally,  a  critical  assessment  of  the  experimental 
results  suggest  an  alternative  design  capable  of  achieving  as  much  as  a  240  times 
improvement  in  the  heat  transport  capacity  in  comparison  to  traditional  micro  heat  pipes. 

Introduction 

Micro  heat  pipe  devices  have  drawn  much  attention  for  their  promise  to  dissipate  high 
power  densities  since  Cotter  (1984)  first  proposed  them.  Their  promise  comes  from  the 
effective  heat  transport  associated  with  phase  change  and  their  small  geometry.  In  general, 
a  micro  heat  pipe  consists  of  a  small  noncircular  channel  that  utilizes  sharp  angled  comer 
regions  as  liquid  arteries  (wick  structure)  to  provide  liquid  flow  from  the  condenser  to  the 
evaporator  region.  The  heat  applied  to  the  evaporator  end  vaporizes  the  working  liquid  in 
that  region,  forcing  it  to  the  cooler  end  where  it  condenses  and  gives  up  the  latent  heat  of 
vaporization.  A  typical  micro  heat  pipe  has  a  hydraulic  diameter  ranging  from  10  pm  to 
several  mm  and  a  length  of  up  to  several  centimeters.  By  definition,  a  micro  heat  pipe 
requires  the  Bond  number  (Bo  =  pgRh  2/cr)  to  be  on  the  order  1  or  less.  Many  studies  have 
been  conducted  on  different  kinds  of  micro  heat  pipe  (Babin  and  Peterson,  1990;  Longtin, 
et  al.,  1994;  Khrustalev  and  Faghri,  1995;  Ma  and  Peterson,  1996,  1998;  Peterson  and  Ma, 
1999).  To  date  the  maximum  heat  flux  dissipated  with  these  devices  has  been  60  W/cm2 
(Cao  el  at.,  1993). 

One  of  the  factors  which  limits  the  heat  transport  capacity  of  micro  heat  pipes  is  the 
large  viscous  losses  associated  with  the  transport  of  the  liquid  from  the  condenser  to  the 
evaporator  (Peterson,  1994  and  Faghri,  1995).  As  the  heat  input  is  increased,  these  losses 
increase.  At  large  heat  inputs,  the  liquid  flow  rate  cannot  sustain  the  evaporation.  When 
this  occurs  the  evaporator  region  is  said  to  ‘dry  out’.  This  limitation  is  called  capillary 
limitation.  Another  shortcoming  is  that  conventional  micro  heat  pipes  are  passive  in 
nature  and  do  not  allow  for  active  control  of  the  system  temperature.  Thus,  as  the  heat 
input  increases,  so  too  does  the  temperature  of  the  heat  source. 

The  present  research  is  premised  on  the  assumption  that  both  augmentation  of  the 
heat  transport  capacity  and  active  thermal  control  of  micro  heat  pipes  can  be  achieved 
through  the  application  of  a  static  electric  field.  When  an  electric  field  is  strategically 


applied  within  micro  heat  pipes  at  the  liquid-vapor  interfaces  that  are  internally  present, 
dielectrophoretic  electrohydrodynamic  forces  can  be  induced  due  to  the  discontinuity  of 
the  electrical  permittivity  of  the  medium.  These  forces  can  contribute  to  a  pressure  jump 
condition  at  the  interface,  effectively  reducing  the  liquid  pressure  in  the  region  affected 
by  the  field.  This  pressure  reduction  can  be  used  to  augment  the  flow  of  liquid  from  the 
condenser  to  the  evaporator  region. 

Jones  (1973)  first  proposed  replacing  the  capillary  wick  structure  in  a  heat  pipe  with 
electrodes.  The  concept  was  demonstrated  successfully  by  Jones  and  Perry  (1974a)  and 
Jones  (1974b).  Their  results  showed  that  the  electric  field  was  able  to  communicate 
working  fluid  between  the  condenser  and  the  evaporator  enabling  the  functioning  of  the 
heat  pipes,  but  the  performance  was  poor  in  comparison  to  the  existing  capillary  driven 
heat  pipes.  Loehrke  and  Debs  (1975)  improved  the  EHD  heat  pipe  of  Jones  and  Perry  and 
were  able  to  achieve  equivalent  thermal  performance  of  conventional  axial-groove  heat 
pipes.  More  recently,  Bryan  and  Seyed-Yagoobi  (1997)  performed  an  experimental  study 
on  monogroove  EHD  heat  pipe.  About  100%  enhancement  in  the  heat  transport  capacity 
was  achieved  in  their  experiments  as  a  result  of  the  EHD  pumping. 

Hallinan  et  al.  (1998)  conducted  experimental  studies  on  electrohydrodynamic 
micro  heat  pipe  arrays.  Their  experiments  were  conducted  on  a  micro  heat  pipe  array 
having  multiple  parallel  1  mm  (wide)  by  0.6  mm  (deep)  grooves.  The  results  of  their 
experiments  showed  that  for  low  heat  input  the  electric  field  had  little  effect  on  heat 
transport  as  capillarity  was  sufficient  to  sustain  liquid  to  the  evaporator.  However,  at  heat 
inputs  sufficient  to  cause  evaporator  dry  out,  the  application  of  an  electric  field  had  the 
effect  of  reducing  the  heat  source  temperature. 

Based  on  the  promising  initial  work  of  Hallinan  et  al.(1998),  more  detailed 
experiments  were  conducted  to  more  completely  identify  the  factors  that  influence  the 
heat  transport  capacity  of  EHD  assisted  micro  heat  pipes.  Additionally,  to  better 
generalize  the  results  for  various  fluids  and  micro  heat  pipe  geometry,  an  analytical 
model  of  an  EHD  assisted  micro  heat  pipe  is  developed  which  can  be  used  to  predict  the 
maximum  heat  transport  capacity  and  the  optimum  geometry. 

EHD  Assisted  Micro  Heat  Pipe  Concept 

Although  the  electrohydrodynamic  heat  pipe  was  introduced  many  years  ago  and 
various  descriptions  have  since  been  studied,  the  application  to  micro  heat  pipe  systems  is 
new.  The  fundamental  operating  principles  of  electrohydrodynamic  (EHD)  micro  heat 
pipes  are  described  as  follows.  When  an  electric  field  is  applied  within  the  heat  pipe,  if 
the  working  fluid  is  a  dielectric,  the  molecules  in  the  working  fluid  are  polarized,  causing 
their  negatively  charged  electron  cloud  to  be  displaced  toward  the  high  voltage  electrode 
relative  to  the  position  of  their  positively  charged  nuclei.  Since  the  kinetic  energy  of  the 
liquid  molecules  is  much  less  than  that  of  vapor  molecules  at  equivalent  temperature,  the 
polarizing  influence  of  the  electric  field  is  substantially  greater  in  the  liquid  phase  than  in 
the  vapor  phase.  Because  all  systems  tend  toward  minimum  energy,  when  an  electric 
field  is  applied  within  heat  pipes,  the  liquid  within  tends  to  fill  the  regions  of  higher 
electric  field  intensity.  If  the  electrodes  are  strategically  placed  in  the  evaporator,  then 
liquid  will  be  drawn  there  while  vapor  will  be  displaced  toward  the  condenser. 

Based  on  this  reasoning,  in  the  present  experimental  model  the  electrodes  were 
located  in  the  evaporator  where  the  working  liquid  is  needed.  These  electrodes  were 


located  along  every  other  groove.  Fig.  1(a)  shows  an  overview  of  the  micro  heat  pipe 
used  in  the  present  model.  As  well,  the  liquid  orientation  within  the  micro  heat  pipe  is 
shown  in  the  absence  of  an  electric  field  with  a  slight  inclination  angle.  Fig.  1(b) 
illustrates  the  liquid  orientation  with  the  electric  field  applied,  with  liquid  filling  the 
grooves  where  the  electrodes  are  present  and  receding  in  the  grooves  without  electrodes. 
For  the  electrode  configuration  employed,  the  localized  electric  field  in  the  evaporator 
draws  liquid  into  the  evaporator.  A  long  wavelength  instability  (Melcher,  and  Smith, 
1969)  causes  the  liquid  to  bridge,  producing  a  "slug"  of  liquid  beneath  the  electrodes.  If 
the  vapor  bubble  trapped  in  the  evaporator  leg  is  small,  liquid  flow  losses  to  the 
evaporator  are  reduced  in  comparison  to  Fig.  l(a)’s  configuration.  As  such,  the  heat 
transport  capacity  of  the  micro  heat  pipe  can  be  increased.  Note  that  the  adjacent  grooves, 
not  covered  by  electrodes,  provide  a  vapor  return  flow  path.  Thus  a  heat  transfer  loop  is 
formed  within  the  EHD  micro  heat  pipe,  potentially  offering  a  very  efficient  means  for 
heat  transport. 

Modeling  of  EHD  Heat  Micro  Pipes 

The  model  developed  here  follows  that  of  Ma  and  Peterson  (1998)  for  a 
conventional  micro  heat  pipe,  differing  primarily  in  the  inclusion  of  dielectrophoretic 
EHD  forces.  The  model  is  subject  to  the  following  assumptions:  steady-state  flow, 
constant  fluid  properties,  a  nearly  horizontal  configuration,  constant  meniscus  radius  of 
curvature  at  a  given  axial  location,  constant  vapor  temperature,  and  constant  surface 
tension.  Also,  only  axial  variations  in  temperature  and  meniscus  shape  are  considered. 

At  steady  state,  the  operation  of  EHD  micro  heat  pipes  requires  that  the  pressure 
potential  provided  by  the  combination  of  electrohydrodynamic  forces  and  capillarity  be 
greater  than  the  sum  of  all  the  pressure  drops  occurring  throughout  the  liquid  and  vapor 
flow  paths.  This  requirement  can  be  expressed  as  (Peterson,  1994) 

APehd  +  APC  >  APV  +  AP[  +  AP  +  +  APsiUg  (1) 

where  AP+  is  the  normal  hydrostatic  pressure  drop  which  can  be  neglected  in  the  present 
model  because  of  the  assumed  horizontal  operation,  and  APsiug  is  the  pressure  gradient  in 
the  liquid  slug  located  beneath  the  electrodes.  Here,  the  electrohydrodynamic  pumping 
pressure  APehd  can  be  found  from  Maxwell  stress  tensor  (Crowley,  1986) 

7S  =eE,Ej (2) 

where  8ij  is  the  Kronecker  delta.  In  the  present  model,  the  main  electric  field  direction  is 
defined  as  y,  so  Ex  and  Ez  are  much  less  than  Ey.  Thus,  the  stress  tensor  is  given  as  Tyy  = 
1/2  E|  E2  =  P.  At  the  interface  in  the  liquid  phase,  P]  =  1/2-ei  E2,  and  in  the  vapor  phase, 
Pv  =  Vi  £v  E2.  Here  the  electric  field  is  assumed  to  be  uniform.  Therefore,  the  pressure 
discontinuity  at  the  liquid-vapor  interface  due  to  the  applied  electric  field  is 

^ EHD  ~&V)E  (3) 

Here  APehd  acts  normally  from  liquid  to  vapor. 

The  pressure  drops  in  the  micro  heat  pipe  array  are  schematically  illustrated  in 
Fig.  2  and  they  are  determined  as  follows.  Based  on  the  Laplace-Young  equation,  APC 
arises  from  the  capillary  pressure  differences  between  the  evaporator  and  condenser  ends 
and  can  be  written  as: 


A  P=  (7 


(4) 


where  rc,e  is  the  minimum  meniscus  radius  of  the  curvature  in  the  evaporator  and  rc>c  is 
the  average  meniscus  radius  of  curvature  occurring  in  the  condenser.  For  the  other  terms 
in  equation  (1),  the  expressions  suggested  by  Ma  and  Peterson  (1998)  are  used.  The 
minimum  meniscus  radius  in  the  evaporator  is  defined  as 
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The  average  meniscus  radius  in  the  condenser  is  set  equal  to  half  of  the  width  of  channel. 
The  vapor  pressure  drop,  APV,  is  given  as  below  based  on  momentum  conservation 
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where  Ci  and  C2  are  coefficients  depending  only  on  the  channel  angle  of  the  groove  and 
the  contact  angle  of  the  liquid,  and  C  may  be  determined  as 
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However,  the  application  of  electric  field  make  the  wick  structure  of  EHD  micro  heat 
pipe  different  from  the  conventional  one,  resulting  in  thicker  film  and  slugs.  Therefore 
the  coefficient  C*  is  determined  based  on  experimental  data. 

As  seen  in  Fig.  1(b)  the  working  liquid  was  drawn  to  the  region  of  electric  field 
forming  liquid  slugs  beneath  the  electrodes.  The  flow  losses  in  this  slug,  APsiug,  can  be 
derived  for  laminar  flow  in  a  square  channel  filled  with  liquid: 

^slug  H 
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where  Lsiug  is  the  length  of  the  slug  on  both  side  of  the  trapped  vapor  bubble. 

Substituting  all  the  above  expressions  into  Eq.  (1)  and  rearranging,  a  relationship 
between  thermal  throughput  q  and  electric  field  intensity  E  is  found  as  follows 
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In  this  equation,  M  can  be  considered  as  a  parameter  depending  upon  the  geometry  of  the 
micro  heat  pipe  and  the  properties  of  the  working  fluid,  and  it  is  expressed  as 
12  PiLslug  fv  ReA  v  pvLeff  flR.eh[plLl0 
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Once  the  applied  electric  field  intensity  and  the  minimum  meniscus  radius  of  the 
curvature  occurring  in  the  evaporator  are  given,  the  maximum  heat  transport  capacity  can 
be  predicted  by  Eq.  (10). 

Based  upon  the  present  EHD  micro  heat  pipe,  the  coefficients  Ci  and  C2  are  equal  to 
1.195  and  0.45  respectively,  which  were  determined  according  to  the  geometry  of  the 
grooves  and  the  contact  angle  of  the  liquid.  C*  in  Eq.  (7)  is  adjusted  to  be  1.68  to  obtain 
best  correspondence  between  the  model  predictions  and  experimental  results.  An  value  of 
16  was  taken  for  fvReh>v  and  15  for  frReh.i  according  to  the  geometry  of  the  duct  (Bejan, 
1995). 

Experiments 

The  micro  heat  pipe  array  used  in  the  present  experiments  was  manufactured  of  a  1 
mm  thick  glass  slide.  It  was  28  mm  long  and  composed  of  7  parallel  grooves  that 
terminated  in  plenums  on  both  ends.  The  grooves  and  plenum  were  machined  via  an 
ultrasonic  milling  process.  The  cross  section  of  each  channel  was  1  mm  wide  by  0.6  mm 
deep  and  the  spacing  between  the  channels  was  1  mm.  Five  such  arrays  were  evaluated  in 
the  experiments. 

Several  shapes  and  locations  of  electrodes  were  utilized.  The  configuration 
producing  the  best  results  is  as  shown  in  Fig.  1.  The  electrode  pair  sandwiches  the  groove 
in  the  evaporator  section  as  shown  in  Fig.  lc.  From  Fig.  lc,  it  can  be  seen  that  the  high 
voltage  electrodes  were  directly  attached  to  the  1  mm  thick  cover  slide  adhered  to  the 
open  side  of  the  grooves;  hence  the  electrode  was  in  intimate  contact  with  the  liquid.  The 
ground  electrodes  were  attached  to  the  exterior  surface  of  the  grooved  slide.  Therefore, 
the  spacing  between  two  electrodes  of  different  levels  is  1  mm. 

Heat  was  input  to  the  heat  pipe  array  via  an  electric  resistance  heater  element  located 
just  at  the  end  of  channels  on  the  evaporator  side  as  shown  in  Fig.  1.  Heat  was  rejected  by 
pumping  5  °C  deionized  water  through  a  cooling  jacket  attached  to  the  micro  heat  pipe 
array  on  the  condenser  side.  The  entire  facility  was  insulated  with  5  cm  of  ceramic  fiber 
insulation  to  ensure  that  all  heat  was  rejected  through  the  cooling  jacket.  Heat  losses  were 
estimated  to  be  less  than  5%  of  the  total  heat  input. 

Five  calibrated  thermocouples  were  attached  to  the  micro  heat  pipe  array  with 
thermally  conducting  epoxy.  These  thermocouples  were  distributed  along  the  central 
channel  with  spacing  of  5.5  mm.  They  were  used  to  monitor  micro  heat  pipe  temperature 
in  the  tests.  Every  test  continued  until  steady  state  was  reached,  usually  taking  24  to  30 
minutes.  The  uncertainty  associated  with  the  temperature  measurement  was  less  than 
±0.2  °C. 

The  working  fluid  used  in  the  experiments  was  pentane,  a  dielectric  fluid  with  high 
vapor  pressure.  The  micro  heat  pipe  array  was  posed  at  9°  tilt  angle,  with  the  evaporator 
elevated  above  the  condenser.  This  initial  bias  was  only  necessary  at  start-up.  The  applied 
electric  field  intensities  were  4,  6,  8  and  9.5  kilovolts  per  millimeter.  Tests  without 
electric  field  were  also  conducted  to  serve  as  a  baseline  for  tests  conducted  with  electric 
field  applied.  The  heat  inputs  considered  were  1  Watt,  1.5,  2,  2.5  and  3  Watts 
respectively  for  each  of  the  five  different  electric  field  intensities.  The  micro  heat  pipes 
were  filled  with  pentane  with  30%,  40%,  52%  and  60%  liquid  volume  fractions  at 
ambient  conditions. 


The  experimental  test  procedure  was  as  follows.  First  the  deionized  water  was 
pumped  through  the  cooling  jacket.  The  electric  field  and  heat  input  were  established 
simultaneously.  The  wall  temperatures  were  monitored  until  steady  state  was  reached,  at 
which  time  the  heat  input  and  electric  field  were  discontinued.  Some  experiments  were 
conducted  with  the  micro  heat  pipe  visible.  The  transparent  structure  permitted 
visualization  of  the  internal  fluid  orientation.  A  Hi-8  camera  was  used  to  provide  a  video 
record  of  the  operation  in  order  to  better  understand  its  functioning. 

Results  and  Discussion 

The  experiments  conducted  were  mainly  focused  on  investigating  the  influence  of 
the  applied  electrostatic  field  on  the  maximum  heat  transport  capacity.  The  experimental 
maximum  heat  transport  capacity  was  determined  in  the  following  way.  For  an  applied 
electric  field  strength  and  heat  input,  the  experiment  continued  until  the  wall  temperature 
reached  a  steady-state  value.  When  the  maximum  heat  transport  capacity  was  reached, 
dry-out  conditions  in  the  evaporator  caused  the  evaporator  temperature  to  rise  rapidly. 
The  model  permitted  prediction  of  the  maximum  heat  transport  capacity  for  various 
electric  field  intensities.  These  predictions  were  determined  using  Eq.  (10).  Both  the 
experimental  results  and  the  model  predictions  for  the  maximum  heat  transport  capacity 
as  a  function  of  the  electric  field  intensity  are  presented  in  Fig.  3.  Good  agreement 
between  predicted  and  experimental  values  is  achieved.  In  Fig.  3  both  experimental  data 
and  model  predictions  show  that  increasing  electric  field  intensity  causes  an  increase  in 
the  maximum  heat  transport  capacity.  Most  significantly  the  experimental  results  suggest 
about  a  sixfold  improvement  in  the  maximum  heat  transport  capacity  at  the  greatest 
electric  field  considered  (9.5  kV/mm),  above  which  corona  discharge  was  evident. 

In  evaluating  the  heat  transport  capacity  using  Eq.  (10),  it  is  found  that  with  an 
increase  in  the  electric  field  intensity  the  proportion  of  the  electrohydrodynamic  potential 
to  the  total  driving  potential  increases.  This  increase  in  heat  transport  capacity  with 
electric  field  is  better  visualized  by  looking  at  the  relative  importance  of  the  EHD  driving 
force  to  capillarity  as  a  function  of  the  electric  field,  as  shown  in  Table  1 .  Here  APtotai 
represents  the  total  driving  pressure.  It  can  be  seen  from  Table  1  that  the  EHD  potential 
becomes  the  dominant  driving  force  at  high  electric  fields,  constituting  87.3  %  of  the 
driving  potential  at  an  electric  field  of  8  kV/mm. 

Having  established  the  accuracy  of  the  model  in  predicting  the  performance  of  one 
EHD  micro  heat  pipe,  the  developed  model  can  also  be  used  to  analyze  geometrical  and 
parametric  variations  of  the  EHD  micro  heat  pipes.  Setting  the  spacing  between  the 
grooves,  d,  to  be  two  times  the  width  of  a  groove,  i.e.  d  =  2w  and  L,  the  length  of  a 
channel,  the  relationship  between  the  maximum  heat  transport  capacity  and  the  groove 
width  is  revealed  to  be  q  w4  and  is  plotted  in  Fig.  4.  In  this  figure  the  electric  field  is  8 
kV/mm  and  the  micro  heat  pipe  length  is  22  mm.  Fig.  4  shows  that  the  maximum  heat 
transfer  capacity  of  EHD  micro  heat  pipes  is  achieved  when  the  groove  width  is 
maximized.  Note  that  this  would  only  be  true  for  Bond  numbers  less  than  unity.  The 
relationship  between  q  and  w  for  zero  electric  field  is  approximately  q  «  w2. 

Interestingly,  EHD  micro  heat  pipes  receive  more  benefit  from  larger  pore  sizes  than 
conventional  micro  heat  pipes. 


The  relationship  between  the  maximum  heat  transfer  capacity  and  the  length  L  is 
also  obtained  at  8  kV/mm  and  is  plotted  in  Fig.  5.  This  figure  shows  the  effect  of  micro 
heat  pipe  length  for  different  groove  widths  on  maximum  heat  transport  capacity. 
Combined  these  results  suggest  that  a  wider  and  shorter  micro  heat  pipe  is  optimum. 
Similar  conclusions  were  drawn  by  Longtin  et  al.  (1994)  for  conventional  micro  heat 
pipes. 

Closer  inspection  of  the  model  predictions  finds  that  the  vapor  pressure  drop  only 
accounts  for  0.42  %  and  the  losses  in  the  liquid  slug  pressure  drop  respectively  the  0.2  % 
of  the  total  pressure.  Nearly  all  of  the  pressure  drops  are  confined  to  the  liquid  film 
surrounding  the  trapped  vapor  bubble  in  the  liquid  supply  channels  (see  Fig  2).  This 
result  indicates  that  if  the  trapped  vapor  bubble  can  be  eliminated,  much  greater  heat 
transport  capacity  can  be  obtained.  This  ideal  condition  is  achieved  if  the  electrodes  fully 
cover  the  liquid  supply  channels  as  shown  in  Fig.  6.  When  electric  field  is  applied,  the 
channel  covered  by  electrodes  is  filled  with  liquid.  In  this  scenario,  the  pressure  drops  are 
as  illustrated  is  Fig.  6b.  In  this  case,  Eq.  (1)  becomes 

APehd  +  APC  >  APV  +  APsiUg  (12) 

The  heat  transport  capacity  is  then  expressed  as 
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A  comparison  of  the  model  predictions  on  heat  transport  capacity  for  a  fully  and  partially 
liquid  filled  evaporator  channel  as  a  function  of  electric  field  is  demonstrated  in  Fig.  7. 
Here,  the  fully  liquid  filled  means  the  whole  channel  beneath  the  electrode  was  fully 
filled  with  liquid  when  electric  field  applied  as  shown  in  Fig.  6;  and  the  partially  liquid 
filled  means  the  channel  beneath  the  electrode  was  partially  (40-50%)  filled  with  liquid 
as  shown  in  Fig.  1(b).  Clearly,  elimination  of  the  entrapped  vapor  bubble  can 
substantially  increase  the  heat  transport  capacity  of  the  micro  heat  pipe.  In  this  case,  up  to 
a  240  fold  increase  in  the  heat  transport  capacity  is  predicted. 

Experimentally,  electrodes  were  fashioned  to  cover  the  whole  liquid  supply 
channels.  The  micro  heat  pipe  was  charged  with  a  50%  liquid  volume  fraction  at  ambient 
conditions  as  shown  in  Fig.6(a).  The  maximum  heat  transport  capacity  measured 
experimentally  was  observed  to  be  substantially  less  than  the  very  optimistic  predictions 
shown  in  Fig.  8.  The  reason  for  the  poor  performance  was  a  function  of  the  manufacture 
of  the  micro  heat  pipes.  As  shown  in  Fig.  1(c)  only  the  high  level  electrode  was  in  direct 
contact  with  the  liquid  in  the  micro  heat  pipe.  A  layer  of  glass  was  present  between  the 
ground  level  electrode  and  the  liquid.  During  operation  many  free  charges  in  the  liquid 
are  advected  toward  the  evaporating  interface.  Those  positive  charges  collecting  on  the 
glass  surface  near  this  interface  are  unable  to  dissipate  through  conduction  through  the 
glass.  Therefore,  an  electric  double  layer  forms  to  effectively  reduce  the  electric  field 
intensity  of  the  evaporating  interface  as  shown  in  Fig.  9,  diminishing  the  driving  force. 
This  phenomena,  however,  is  not  indicative  of  the  inherent  failure  of  this  concept.  Rather 
it  suggests  a  remedy  that  should  be  addressed  in  future  design.  The  high  voltage  and 


ground  level  electrodes  in  the  evaporating  channels  should  be  in  intimate  contact  with  the 
liquid.  As  well,  the  vapor  return  line  should  be  completely  shielded  from  the  electric 
field.  Such  a  design  is  expected  to  be  the  optimum  one. 

Conclusions 

A  mathematical  model  to  predict  the  maximum  heat  transport  capacity  of 
electrohydrodynamic  micro  heat  pipes  versus  electric  field  intensity  is  presented.  Good 
agreement  is  achieved  in  comparing  the  predicted  value  with  the  experimental  data. 
Equation  (10)  indicates  that  increasing  applied  electric  field  can  increase  the  maximum 
heat  transport  capacity  of  EHD  micro  heat  pipes  since  at  high  level  electric  field,  the 
EHD  driving  force  became  dominant.  The  optimum  geometry  of  the  EHD  micro  heat 
pipe  can  also  be  predicted  using  Eq.  (10).  The  prediction  indicates  that  the  optimum  EHD 
micro  heat  pipe  is  short  and  wide  with  the  requirement  that  the  Bo  number  be  less  than 
one.  An  ideal  configuration  of  the  EHD  micro  heat  pipe  is  predicted  by  Eq.  (10).  The 
ideal  EHD  micro  heat  pipe  model  should  have  specially  shielded  vapor  return  channels 
for  the  vapor  to  reach  condenser  rapidly.  Eq.  (10)  can  be  used  to  predict  the  performance 
of  EHD  micro  heat  pipes  and  aid  the  future  design  of  the  pipes. 


Nomenclature 

A  =  cross-sectional  area,  m2 
a  =  channel  depth,  m 
C  =  coefficient 

d  =  width  of  pipe  +  space  between  pipes,  m 
E  =  electric  field  intensity,  V/m 
f  =  friction  factor 

g  =  gravitational  acceleration,  m/s2 

L  =  length,  m 

U,T  =  effective  length,  m 

M  =  parameter  defined  in  Eq.  (10) 

P  =  pressure,  N/m2 

Q  =  total  heat  input  into  device,  Watts 
q  =  heat  transport  for  one  channel.  Watts 
Re  =  Reynolds  number 
T  =  temperature,  K 
w  =  channel  width,  m 
£  =  permittivity 
p  =  density,  kg/m3 
a  =  surface  tension,  N/m 
X  =  latent  heat  of  vaporization,  J/kg 
|i  =  viscosity,  Ns/  m2 


V|/  =  tilt  angle 


Subscripts 

c  =  capillary 

EHD  =  electrohydrodynamic 

e  -  evaporator 

h  =  hydraulic 
l  =  liquid 
max  =  maximum 
v  =  vapor 
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Figure  1 .  The  schematic  configuration  of  the  EHD  micro-heat  pipes 

(a)  The  configuration  without  application  of  electric  field 

(b)  The  configuration  with  application  of  electric  field 

(c)  Schematic  of  the  sandwich  electrode  structure  at  A-A  section 
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Figure  2.  Elustration  of  the  pressure  drops  distribution 
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Figure  3  The  predicted  and  the  experimental  maximum  heat  transport  capacity 
versus  electric  field  intensity 


Figure  4  Theoretical  maximum  heat  transport  capacity  for  a  single  channel  versus 
groove  width  with  and  without  an  applied  electric  field  of  8  kV/mm 


Figure  5  Theoretical  maximum  heat  transport  capacity  versus  length  of  the  pipe 
with  various  groove  widths  at  8  kV/mm 
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Figure  6  (a)  Liquid  fully  filled  schematic  (the  ideal  model) 
(b)  Pressure  drops  in  liquid  fully  filled  case 
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Figure  7  Predictions  of  maximum  heat  transfer  capacity  (log  value)  versus  electric 
field  for  a  fully-liquid  filled  and  partially  liquid  filled  evaporator 


Figure  8  The  maximum  heat  transport  capacity  (log  value)  for  fully-liquid 
filled  evaporator  channel,  prediction  and  the  experiment  data 
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Table  1 .  The  ratio  of  EHD  driving  potential  and  capillary  pressure  potential  to  the 
total  driving  potential 


Electric  field  intensity  (kV/mm) 

AP electr  ! AP total 

AP Cap  / AP total 

0 

0 

100% 

4 

63.3  % 

36.7  % 

6 

79.5  % 

20.5  % 

8 

87.3  % 

12.7  % 

9.5 

90.7  % 

9.3  % 

Active  Thermal  Control  of  An  Ion-drag  Pump  Assisted  Micro 
Heat  Pipe 


Abstract 

An  ion-drag  pump  is  utilized  to  enhance  the  heat  transport  capacity  of  micro  heat 
pipes.  An  analytical  model  is  developed  to  estimate  the  maximum  heat  transport  capacity 
as  a  function  of  the  applied  electric  field.  The  model  demonstrates  that  four  times 
increase  in  heat  transport  capacity  is  made  possible  through  application  of  an  electric 
field.  The  model  predictions  are  shown  to  match  experimental  data  well.  Given  the 
potential  for  augmentation  of  heat  transport  capability,  a  transient  analytical  model  was 
developed  to  permit  variation  of  the  electric  field  with  applied  thermal  load.  A 
proportional-integral-derivative  controller  was  used  to  simulate  active  temperature 
control.  The  feasibility  of  achieving  the  active  temperature  control  was  demonstrated 
experimentally. 


Introduction 

Effective  thermal  management  for  high  temperature  electronics  are  becoming  more 
concerned.  High  temperature  can  cause  electronics  performance  to  degrade  or  cause  it  to 
fail  prematurely.  Faster  speeds  generate  more  heat,  while  smaller  packages  increase  the 
difficulty  of  removing  it  from  the  devices  that  are  adversely  affected  by  increase 
temperature.  The  surface  heat  flux  at  the  chip  level  for  a  5  mm  by  5  mm  chip  dissipating 
10  W  is  4  x  10 5  W/m2  (McCluskey  et  al.  1997).  Developing  a  high  intensity  heat 
dissipating  system  in  micro  scale  is  of  significance  for  meeting  not  only  the  performance 
requirements,  but  reliability  requirements  as  well.  Micro  heat  pipe  has  shown  its  distinct 
feature  in  the  passive  management  of  the  heat  dissipation.  However,  the  requirement  of 
extremely  high  power  dissipating  of  the  high  speed  running  electronics  and  the 
requirement  of  dynamic  control  of  the  temperature  make  the  conventional  micro  heat 
pipe  incompetent.  Electrohydrodynamically  assisted  micro  heat  pipe  becomes  noticeable 
due  to  its  high  capability  of  intensity  heat  dissipation  and  active  controllable  nature. 
Among  electrohydrodynamic  pumps,  ion-drag  pump  shows  increasing  potential.  The 
configuration  of  ion  drag  pump  is  illustrated  in  Fig.  1  and  the  operation  of  the  pump  is  as 
following.  As  Fig.  1  shown,  an  electric  field  is  established  between  the  negatively 
energized  emitter  and  the  grounded  collector.  Free  charges  are  generated  within  the 
dielectric  fluid  by  emission  of  the  emitter  as  a  corona  source  and  then  are  dragged 
towards  the  collector.  The  friction  between  the  moving  charges  and  the  working  fluid 
creates  a  driving  potential  thus  setting  the  fluid  in  motion. 

Most  of  theoretical  studies  of  the  ion  drag  pump  are  based  on  the  analyses  provided 
by  Stuetzer  (1959).  Stuetzer  established  a  theoretical  model  based  on  Gauss’  Law  and 
derived  the  relation  between  pressure  and  electric  field  intensity 

P-P0  =  (e/2)-(E2  -  E02) 


(1) 


where  P  is  the  pressure  within  fluid,  8  the  permittivity  and  E  the  electric  field  intensity. 
Based  upon  this  relation  he  also  developed  the  relations  between  operating  voltage  and 
pressure  for  cylindrical  and  spherical  coordinate  systems.  Experiments  were  conducted  to 
measure  the  pressure  differences  using  corona  point  emitter  as  an  ion  source.  It  was 
found  that  for  a  corona  point  ion  emitter  the  largest  pressure  buildup  takes  place  in  the 
immediate  environment  of  the  point.  This  model  approximately  matched  the  results  of  the 
experiments  for  static  conditions.  Stuetzer  (1960)  later  discussed  the  dynamic  situation. 
He  used  a  loop  system  to  measure  the  pressure  drops  within  the  moving  fluid  and  derived 
an  approximate  expression  relating  dynamic  pressure  to  the  constant  operating  voltage 
involving  pressure  drop.  The  power  efficiency  of  the  pump  system  was  also  discussed. 

Later,  Pickard  (1963a)  provided  more  insight  on  the  efficiency  of  pressure 
generation  of  ion  drag  pump  and  proposed  the  concept  of  “cut-in”  voltage.  He  considered 
hydrodynamic  instabilities  as  well.  He  concluded  that  below  the  “cut-in”  voltage  the 
pump  is  inoperative  and  above  it  the  pump  generated  a  pressure  of  the  form  P~(V-V0)2. 
He  demonstrated  that  pressure  instabilities  occurred  for  constant  voltage  due  to  random 
fluctuations  in  emission. 

Pickard  (1963b)  performed  experiments  to  investigate  the  factors  influencing  ion 
pumping.  It  was  found  that  suitable  electrode-liquid  combinations  could  produce  better 
results.  Acetone  turned  to  be  suitable  liquid  because  it  possesses  a  relatively  high 
permittivity  (8  =  20)  and  it  can  be  readily  purified.  Molybdenum  was  found  to  be  a 
suitable  electrode  material.  The  pressure  instability  occurring  in  his  experiments  was 
found  to  lead  to  undesirable  short-term  oscillations.  He  suggested  three  causes  for  the 
instability:  the  formation  of  a  vapor  pocket,  hydrodynamic  effects  of  the  fluid  velocity¬ 
generated  pressure  interrelationship  and  variations  in  pumping  efficiency.  A  damping 
plug  was  inserted  to  the  pump  channel  to  damp-out  the  instability.  Ac  pumping  was  also 
tested,  and  he  concluded  that  pumping  was  normally  a  multiionic  process  and  that  to 
obtain  a  maximum  pressure  head  it  was  necessary  to  insure  one  electrode  generates  more 
charges  than  the  other.  The  pitted  anode  surface  in  his  experiments  indicated  a  strong 
localized  emission  of  ions  in  the  experiments.  Therefore,  he  discussed  the  electrode  aging 
effects.  In  addition,  the  results  of  his  experiments  demonstrated  that  the  measured  current 
and  the  magnitude  of  the  observed  pumping  bore  no  necessary  connection.  This  is 
different  from  the  Stuetzer’ s  conclusion. 

The  basic  assumption  in  the  derivation  of  Stuetzer’ s  and  Pickard’s  model  is  using 
two  plane-parallel  electrodes  to  approximate  the  wire  grid  electrodes  in  practice.  This 
means  the  variation  of  charge  density  around  point-shaped  emitter  was  neglected.  Finally, 
they  derived  the  maximum  pressure  as 
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where  d  is  the  spacing  of  the  electrodes. 

Based  on  Gauss’s  Law,  Woodson  and  Melcher  (1968)  derived  an  expression 
describing  charge  density  variation  with  the  distance  of  emitter  and  electric  Reynolds 
number.  Using  the  Maxwell  stress  tensor,  the  pressure  difference  between  the  electrodes 
was  given  by  P(d)  -  P(0)  =  l/2eo  [E2(d)-  E2(0)],  and  finally  yielded 
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Here  Re  =  £o  vo/(p0bd)  is  the  electric  Reynolds  number,  where  v0  is  the  velocity  of 
charges  at  the  emitter  and  po  is  the  initial  charge  density  at  the  emitter.  They 
demonstrated  that  the  electric  pressure  was  limited  by  the  breakdown  strength  of  the 
dielectric  medium.  Hence  they  suggested  that  magnetic  field  was  able  to  provide  higher 
pressure  than  electric  field. 

Crowley  (1986)  introduced  a  basic  understanding  of  ion-drag  force  from  an 
electrostatics  point  of  view.  In  light  of  the  Coulomb’s  force  present,  he  developed  a  force 
density  inside  ion-drag  pump: 


F  = 


PqVq  | 
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and  the  net  force  developed 
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Then  he  derived  the  maximum  pressure  as 


P  = 


2eVn 


(N/m2) 


(5) 


It  is  noted  that  the  maximum  pressure  given  by  Eq.(5)  is  much  larger  than  that  of  Stuetzer 
and  Packard  given  by  Eq.(2).  Eq.  (5)  is  also  an  approximation  since  many  assumptions 
were  made.  And  the  primary  assumption  is  that  the  electric  field  is  uniform.  In  practice, 
the  electric  field  is  distorted,  which  makes  problems  more  complicated.  Eq.  (4)  indicates 
that  a  nonuniform  pressure  exists  within  fluid  because  of  the  variation  of  force  density  F 
with  the  distance  x  from  the  emitter.  Therefore,  an  internal  stirring  of  the  fluid  exists, 
forming  electroconvection,  which  can  reduce  the  output  pressure. 

Babin,  Peterson  and  Seyed-Yagoobi  (1993)  developed  an  analytical  model  by 
simply  applying  Crowley’s  formula  APehd  =  PoV  to  the  expression  APC  +  APehd  ^  APV  + 
APi  +  APg  to  evaluate  the  heat  transport  capacity  of  an  ion  drag  assisted  micro  heat  pipe. 
Here  APC  is  the  capillary  pumping  pressure,  APehd  is  the  pumping  pressure  provided  by 
the  ion-drag  pump,  APV  is  the  pressure  required  to  transport  the  vapor  from  the  evaporator 
to  the  condenser,  APi  is  the  pressure  required  to  transport  the  liquid  from  the  condenser  to 
evaporator  and  APg  is  the  hydrostatic  pressure  due  to  the  local  gravitational  acceleration. 
Charge  density  po  is  determined  by  po  <  (2 eV/d2),  which  is  Crowley’s  formula  described 
in  Eq.  (5).  Only  the  maximum  pumping  pressure  can  approximately  be  estimated  through 
this  relation.  They  constructed  a  two-stage  ion-pump  and  conducted  an  experimental 
investigation  to  verify  their  analytical  model.  Both  the  results  predicted  by  the  model  and 
obtained  from  experiments  showed  that  the  pressure  generated  by  the  ion-drag  pump 
increased  linearly  with  increasing  voltage.  Their  results  also  showed  that  the  slope  of  the 
predicted  curve  is  approximately  10%  less  than  that  of  the  experimental  curve.  However, 
they  did  not  establish  the  relation  between  heat  transport  capacity  and  operating  voltage. 

More  recently,  Bryan  and  Seyed-Yagoobi  (1994)  built  an  ion-drag  pump  with 
screen  type  electrodes  occupying  entire  pipe  cross  section.  Although  they  use  Maxwell 
shear  stress  to  derive  two-dimensional  electrical  force,  the  final  developed  model  is  one 
dimensional.  Combining  the  derived  electric  force  with  the  momentum  equation  and 


Gauss’s  law,  after  non-dimensionalizing,  they  solved  the  equations  using  a  semi-implicit 
method.  The  profile  of  bulk  liquid  velocity  via  applied  voltage  was  obtained.  Three 
nondimensional  electrical  parameters  were  used  in  their  analyses.  They  are  electric 
Reynolds  number  ReE  =  eU/(ceL),  electric  source  number  Esi  =  bVo/(LU)  and 
Es  =  poL2/£Vo.  Here  L  is  the  spacing  of  electrodes.  The  electric  field  distribution  along 
the  pump  axis  was  determined.  It  was  found  that  electric  source,  Es,  considerably 
influences  the  electric  field  distribution  since  it  represents  the  magnitude  of  the  charge 
injection  at  the  emitter.  As  Es  increases  the  gradient  of  the  electric  field  increases, 
resulting  in  an  increased  flow  rate. 

Later,  Seyed-Yagoobi,  el  al.  (1995)  further  refined  their  analytical  model.  Their 
results  were  presented  in  terms  of  three  nondimensional  electric  parameters,  ReE,  Esi  and 
Es.  The  results  demonstrated  that  higher  electric  Reynolds  number,  Res,  made  a  more 
uniform  charge  density  distribution  due  to  the  longer  charge  relaxation  time  te  =  £/ae 
(recall  ReE  =  eU/(aeL)).  However,  the  data  showed  that  electric  slip  number  Esi  played 
the  most  critical  role  in  the  distribution  of  charge  density.  Higher  Ren  also  resulted  in 
larger  gradient  of  nondimensional  electric  field.  All  the  results  showed  low  electric  slip 
number  made  efficient  operation  of  an 

ion-drag  pump.  Further  insight  indicated  that  low  conductivity  and  mobility  and  high 
permittivity  resulted  in  high  efficiency. 

The  present  research  applies  the  results  of  previous  researcher  to  the  development 
of  a  theoretical  model  capable  of  predicting  the  performance  and  the  operational 
characteristics  of  an  ion  drag  pump  enhanced  micro  heat  pipe  device.  An  active  feedback 
control  model  is  also  developed  to  simulate  the  use  of  the  ion  drag  pump  as  a  vehicle  for 
achieving  active  thermal  control  of  heated  surfaces  subject  to  transient  thermal  loads. 


Theoretical  Prediction  of  Maximum  Heat  Transport 

The  assumptions  employed  in  the  model  development  are:  steady-state  flow,  one  - 
dimensional  electric  field  distribution,  only  one  kind  of  ion  (positive),  constant  fluid 
properties,  constant  surface  tension  and  perfectly  wetting  liquid.  Turbulence  in  the  fluid 
is  neglected. 

The  geometry  of  the  ion-drag  pump  assisted  micro  heat  pipes  is  shown  in  Fig.  1 .  In 
Fig.  1  (a),  the  collector  is  connected  to  ground  level  and  the  emitter  to  the  negative  high 
voltage  level.  When  the  electric  field  is  applied,  ions  are  generated  in  the  neighborhood 
of  the  emitter  within  the  dielectric  liquid  by  emission  of  a  discharge  corona.  The  collector 
drags  the  ions  and  induces  the  motion  of  the  liquid  through  the  friction  between  moving 
ions  and  liquid  forming  the  ion-drag  pump.  As  shown  in  Fig.l  (b),  the  micro  heat  pipes 
are  grooved  in  radial  direction  on  the  top  surface  and  covered  by  a  circular  heated  disc. 
The  channels  have  a  square  cross  section.  The  vertical  supply  channel  in  which  the  ion- 
drag  pump  is  located  possesses  circular  cross  section.  The  ion-drag  pump  pushes  the 
working  liquid  through  the  central  hole  into  the  micro  heat  pipes,  as  shown  in  Fig.  1  (c). 
Heat  is  added  on  the  top  of  the  device  evaporating  the  working  fluid.  Condensation  takes 
place  in  the  surrounding  annular  region. 

For  the  horizontal  section  of  the  device,  i.e.  micro  heat  pipe,  the  liquid  control 
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volume  is  selected  as  shown  in  Fig.  2.  Here,  m  is  the  mass  flow  rate  and 
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The  conservation  of  momentum  in  x  direction  for  the  liquid  control  volume  requires 
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where  m  is  the  momentum  flux  due  to  phase  change  at  the  interface.  The  frictional  force 
in  the  horizontal  channels  is  determined  using  the  geometry  in  Fig.  3.  The  perimeter  of 
the  liquid  phase  is  8r,  noting  that  r  is  the  function  of  x.  Here  x  =  0  is  at  the  central  point 
of  the  top  evaporator  region  and  x  =  Lh  is  at  the  condenser  end  of  micro  heat  pipe.  After 
some  manipulation  and  neglecting  higher  order  differential,  Eq.  (7)  becomes 
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dx  dx 

Note  that  the  first  term  includes  the  capillary  pressure,  expressed  as 
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In  the  present  model,  the  primary  portion  of  the  pressure  P  results  from  the 
electrohydrodynamic  pumping  force.  Therefore,  the  derivation  of  the  relation  between 
liquid  driving  pressure  and  EHD  pumping  force  is  of  primary  significance.  The 
configuration  of  control  volume  in  vertical  section  of  the  pump  is  illustrated  in  Fig.  4a. 
Here  Fe  is  the  electrohydrodynamic  pumping  force,  Fw  is  the  gravitative  force  and  Ff  is 
the  frictional  force.  Assuming  that  there  are  N  micro  grooves  and  considering  the  vertical 
section  and  horizontal  section  as  an  whole  control  volume,  a  force  balance  is  established 
on  the  control  volume.  The  pumping  pressure  within  the  liquid  flow  is  thus  expressed  as 
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Here  Apis  the  vertical  channel  cross-section  area,  7tDp2/4.  The  gravitative  and  frictional 
force  are  respectively  defined  as 
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Ff  =  Tw-7tDpLp  (12) 

The  electrohydrodynamic  pumping  force,  Fe,  is  derived  from  the  Maxwell  stress 
tensor.  It  should  be  noted  that  the  force  acting  on  the  charges  is  not  just  the  form  as  might 
be  expected  from  the  basic  form  of  electrostatic  force  equation:  Fe  =  pE.  This  equation  is 


valid  only  for  a  vanishingly  small  volume  of  charges  in  an  externally  generated  field.  In 
present  model,  the  charge  distribution  is  not  clear.  Therefore,  Maxwell  stress  tensor  is 
used,  which  is  valid  for  any  charge  distribution. 


Tij  =  eEiEj 


(13) 


where  5jj  is  the  Kronecker  delta.  The  main  electric  field  direction  is  in  the  vertical(y) 
direction,  so  the  magnitude  of  Ex  and  Ez  are  much  less  than  that  of  Ey.  Thus  the  stress 
tensor  is  given  as  Tyy  =1/2  ei  E2  =  P.  Therefore,  Fe  is  defined  as 
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Substituting  Eq.  (1 1),  (12)  and  (14)  into  Eq.  (10)  and  rearranging  yields 
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To  compute  the  friction  shear  stress  xw ,  the  liquid  flows  in  both  the  vertical  section  and 
horizontal  section  are  assumed  to  be  similar  to  fully  developed  duct  flow.  Thus  the  shear 
stress  is  expressed  as  (Bejan,  1995) 

*w=^P,U?f  (16) 

where /=  k/Re.  Here  k  is  the  parameter  depending  on  the  geometry  of  the  channel  (see 
Bejan,  1995).  Reynolds  number  Re  =  piUitPDp/jj.i  •  Using  the  relation  of  the  mass  flow 

rate  in  the  vertical  channel  m  p  -  piUi,p  Ap,  Eq.  (16)  becomes 
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For  the  horizontal  section  of  the  pump,  the  cross-sectional  area  of  liquid  flow  can  be 
determined  in  the  light  of  Fig.  3 


Ah  =  (4-tt)-r 

The  hydraulic  radius  rh,h  is  defined  as 
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where  ph  is  the  wetted  perimeter  of  the  liquid  flow. 

Analogous  to  the  derivation  of  Eq.  (17)  and  using  Eq.  (19),  noting  that 
Re  =  pi  U|,p  Dh.h  /(Xi ,  where  D^h  is  the  hydraulic  diameter  of  the  liquid  flow  and 
fh.h  =  Dh.h  /2,  the  frictional  shear  stress,  Tw,  for  the  horizontal  section  is  expressed  as 
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Substituting  Eq.  (15),  (17),  (18)  and  (20)  into  (8)  and  using  the  relation  of  the  mass  flow 


rate  in  horizontal  channel  m  h  =  pt  Uh  Ah  and  m  =  q”/  hfg,  where  q  is  the  heat  flux, 


yields  the  following  equation 
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In  order  to  obtain  nondimensional  form,  following  nondimensionalizations  are  employed: 
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Substituting  Eq.  (22)  into  (21)  yields 
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In  order  to  solve  Eq.  (23)  numerically,  the  following  assumptions  and  boundary 
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conditions  are  employed:  (1 )  m  h  =  m  p  /N  at  x  =  0;  (2)  m  p  is  initially  estimated  by 
considering  a  force  balance  on  the  vertical  channel. 
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(3)  At  the  maximum  heat  transport  condition  dry-out  occurs  at  x  =  L,  and  thus  r  =  w/2  at 
x*  =  0  and  r*  =  0  at  x*  =  1 ;  (4)  Pv  is  determined  as  the  saturate  pressure  associated  with  the 
local  temperature.  Finally  Eq.  (23)  is  solved  subject  to  a  conservation  of  mass  constraint, 
requiring 
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A  solution  to  this  system  of  Equations  is  obtained  using  fourth  order  Runge-Kutta 


method.  Newtons’  method  is  employed  to  iterate  on  m  p  convergence  is  achieved. 


Experiment 

The  schematic  experimental  apparatus  is  shown  in  Fig.  1 .  The  diameter  of  the 
vertical  is  0.003175  m.  the  spacing  of  the  electrodes  is  0.01  m.  The  width  of  the  micro 
heat  pipe  is  0.001  m  and  the  length  is  0.005  m.  There  were  eight  such  channels.  The 
working  fluid  used  is  n-Pentane.  The  amount  of  the  working  liquid  is  large  enough  to 
maintain  its  temperature  constant  in  the  experimental  process,  that  is,  room  temperature. 
The  temperature  of  the  heat  source  was  measured  using  thermistor.  Heat  was  input  via  an 
electric  resistance  heater  element  located  on  the  top  of  the  system.  The  applied  electric 
field  was  fixed  at  15  kV/mm  level. 

Results  and  Discussion 

The  maximum  heat  flux  of  the  experiment  at  15  kV/mm  is  40  W/cm2  at  steady  state, 
which  is  obtained  for  the  maximum  electric  field  achievable  before  electric  breakdown  is 
shown. 

Sample  model  calculations  are  presented  in  Figs.  5  and  6.  It  should  be  recognized 

that  for  a  given  electric  field  E,  the  maximum  heat  input  flux  q  is  achieved  when  both  r 
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and  m  reach  zero  at  x*  =  1 ,  since  under  these  conditions  all  working  liquid  within  the 
micro  heat  pipes  is  evaporated  into  vapor  without  dryout  occurring.  The  predicted 
maximum  q  is  presented  as  a  function  of  electric  field,  as  shown  in  Fig.  7.  This  figure 
clearly  indicates  that  the  applied  electric  field  plays  a  significant  role  in  the  performance 
of  ion-drag  assisted  micro  heat  pipe.  Increasing  electric  field  intensity  can  increase  the 
maximum  heat  transport  capability  noticeably.  A  similar  trend  was  observed  by  Bryan 
and  Yagoobi  (1994).  Experimental  results  is  also  depicted  in  Fig.  7.  Good 
correspondence  prediction  and  experiment  is  obtained 

The  geometry  of  a  micro  heat  pipe  generally  has  significant  influence  on  its  heat 
transport  capability.  Therefore,  predictions  are  made  to  demonstrate  the  influence  of 
geometry  variations  of  the  heat  transport  capacity.  These  results  are  presented  in  Figs.  8 
and  9.  Fig.  8  shows  that  the  width  of  micro  heat  pipe  channels  has  a  large  influence  on 
the  maximum  heat  transfer  capability.  It  is  suggested  that  roughly  q"max 06  w3,  which 
agrees  with  the  results  from  other  studies  conducted  by  Longtin  et  al.  (1994)  and  Yu  et  al. 
(2000).  Fig.  9  illustrates  the  effect  of  the  length  of  micro  heat  pipe.  It  is  seen  that  the 
maximum  heat  transfer  capacity  decreases  linearly  with  an  increased  micro  heat  pipe 
length.  These  results  further  confirm  the  recommendations  of  Longtin  et  al  that  short  and 
wide  channel  produces  optimum  heat  transfer. 

The  effect  of  the  spacing  of  the  electrodes  and  the  diameter  of  the  vertical  channels 
on  the  liquid  driving  pressure  of  the  micro  heat  pipe  are  evaluated.  The  results  do  not 
show  noticeable  influences  from  them.  Therefore,  these  effects  are  negligible. 

Controlling  of  the  EHD  Pump 

The  most  promising  feature  of  the  EHD  pump  to  the  micro  heat  pipe  is  the 
application  of  active  temperature  control.  The  feasibility  of  heat  transport  capability 
augmentation  demonstrated  by  previous  results  leads  the  possibility  of  achieving  thermal 
control  for  transient  heat-loadings.  The  control  is  implemented  through  the  control  of 
electric  field  intensity  and  thereby  controlling  the  mass  flow  rate  of  the  working  liquid. 
The  governing  equation  is 
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A 


(lit) input  -q(t)evap  =  Cp 

where  qinput  is  the  heat  input,  qevap  is  the  heat  to  be  evaporated,  and  mb  and  Cp  are 
respectively  the  mass  and  specific  heat  of  the  heated  substrate.  For  evaporation,  qevap  can 
be  expressed  as 


<leVap=mhfg  (27) 

In  the  EHD  micro  heat  pipe,  the  mass  flow  rate  m  is  determined  by  the  applied  electric 


field  E,  namely,  m  -m  (E).  Thus  Eq.  (26)  can  be 

input  -m(E)hfg  =Cp^(mbT) 
Applying  the  Laplace  transformation  on  Eq.  (28)  yields 

q(s)  -  m(E)  ■  h  fg 


T(s)  =  - 


Cpmhs 


(28) 


(29) 


where  s  is  the  Laplace  variable.  This  latter  relationship  presumes  that  the  flow  field 


responds  much  more  quickly  than  the  temperature  filed  since  it  uses  m  ( E )  as  calculated 
using  Eq.  (23).  Eq.  (29)  is  solved  using  the  Matlab  simulink  environment  and  the  control 
block  diagram  is  presented  in  Fig.  12.  The  control  operation  procedure  is  as  following. 
The  temperature  of  the  heat  source  is  measured  and  compared  with  the  reference 
temperature.  If  the  temperature  of  the  heat  source  is  higher  than  the  reference  one,  the 
micro-controller  will  turn  on  the  EHD  pump  at  certain  electric  field  level.  Certain  amount 
of  working  fluid  will  be  provided  into  the  evaporator  region  to  cool  down  the  region  and 
obtain  the  desired  temperature  in  the  region.  Once  the  heat  source  temperature  is  lower 
than  the  reference  one,  the  micro-controller  will  turn  off  the  EHD  pump. 


Simulation  and  Experimental  Results  with  Active  Thermal  Control 

The  corresponding  transient  response  is  presented  in  Fig.  13.  It  is  seen  that  after 
several  oscillations  (120  seconds),  the  temperature  reached  steady  state  at  required 
temperature.  The  heat  input  was  chosen  as  a  pulse  pattern  in  order  to  simulate  the  actual 
situation  in  electronic  components.  The  duty  cycle  was  selected  from  70%  to  100%.  The 
heat  input  control  was  simulated  using  Pspice  software,  and  the  diagram  of  the 
schematics  was  shown  in  Fig.  14.  After  all  the  control  models  were  verified  by  computer 
simulation,  the  physical  facility  was  built  up  and  the  experiment  was  implemented. 

Experiments  were  performed  on  the  controlled  ion-drag  pump.  The  electric  field 
was  generated  by  a  high  voltage  power  supply  (Model:  PS/LG- 10  P  15).  The  heat  was 
generated  by  a  resistance  heater  and  controlled  through  the  control  circuit.  The  reference 
temperatures  were  set  to  be  63  °C,  68  °C,  87  °C  and  89  °C  respectively.  The  P.  I.  D. 
controller  controls  the  on-off  duty  cycle.  When  the  temperature  reached  above  the 
reference  temperature,  the  micro  controller  would  turn  on  the  electric  field,  and  ion  drag 
pump  was  started  to  cool  down  the  heat  source. 


Two  of  typical  plots  that  recorded  the  variations  of  temperature  versus  time  under 
control  condition  are  presented  in  Figs.  13  and  14.  It  can  be  seen  that  the  heat  control 
system  effectively  controlled  the  temperature  of  the  system  according  to  the  reference 
temperature.  Fig.  14  shows  that  if  no  electric  field,  the  temperature  will  go  up  fast  and 
reaches  very  high  level.  While  the  applied  electric  field  makes  the  temperature  go  down 
quickly  and  keep  it  at  lower  level.  It  is  also  noted  that  although  the  heat  input  (26  W/cm2) 
in  Fig.  13  is  higher  than  that  in  Fig.  14  (17  W/cm2),  the  control  system  can  control  the 
temperature  depending  on  the  requirement.  The  control  of  the  magnitude  of  the 
temperature  oscillation  can  be  achieved  by  adjusting  the  duty  cycles  of  the  electric  field 
pulse,  that  is,  the  ratio  of  on  time  to  off  time.  It  can  be  seen  from  Fig.  14  that  through 
enlarging  the  duty  cycle  or  increasing  the  frequency  of  the  duty  cycle,  the  electric  field  is 
applied  intensively.  The  temperature  is  thereby  almost  kept  as  constant,  90.5  ±  0.5  °C. 
Therefore,  the  control  of  the  temperature  can  be  carried  out  more  precisely.  In  the 
comparison  of  Fig.  13  with  Fig.  14,  it  is  apparent  that  high  frequency  can  make  the 
control  of  temperature  more  precise  than  the  effect  of  long  duty  cycle.  Hence  high 
frequency  is  demanded  in  precise  case. 


Conclusion 

An  analytical  model  is  developed  on  the  ion  drag  pump  assisted  micro  heat  pipe. 
The  interfacial  effects,  the  interfacial  radius  of  curvature  and  the  change  of  flow  direction 
are  taken  into  account  in  the  model  development.  Nondimensionalization  are  employed. 
The  resulting  differential  equation  is  solved  numerically.  The  model  predictions  show 
that  the  heat  transport  capacity  increases  with  the  increase  of  the  intensity  of  the  electric 
field.  This  is  due  to  the  increase  of  both  the  amount  of  emitted  free  charges  and  the 
velocity  of  these  charges.  The  evaluation  of  the  ion  drag  pump  geometry  does  not  show 
significant  influence  on  the  performance.  Some  of  the  model  predictions  agree  with  the 
experimental  results  from  other  studies. 

Temperature  control  is  achieved  by  the  control  of  electric  field.  This  is  the 
noticeable  feature  of  this  system  that  distinct  from  the  other  heat  dissipating  system.  A 
desired  temperature  is  maintained  through  controlling  the  electric  field.  High  frequency 
can  make  a  precise  temperature  control. 


Nomenclature 

A  =  cross-sectional  area,  (m2) 
b  =  the  mobility  of  the  ion 
Cp  =  specific  heat  of  the  fluid,  (J/kg-K) 
D  =  diameter  of  the  channel,  m 
d  =  the  spacing  of  the  electrodes,  m 
E  =  electric  field  intensity,  (V/m) 
f  =  friction  factor 

g  =  gravitational  acceleration,  (m/s2) 
hfg  =  latent  heat  of  vaporization,  (kJ/kg) 
L  =  length,  m 


m  =  mass  of  the  working  fluid,  kg 

• 

m  =  mass  flow  rate,  (kg/s) 

P  =  pressure,  (N/m2) 
q  =  heat  transport,  Watts 
Re  =  Reynolds  number 

r  =  the  curvature  radius  of  liquid-vapor  interface,  m 

T  =  temperature,  K 

w  =  channel  width,  m 

U  =  the  velocity  of  flow,  m/s 

V  =  the  operating  voltage,  Volt 

e  =  permittivity 

p  =  density,  (kg/m3) 

a  =  surface  tension,  N/m 

ce  =  the  conductivity 

X  =  latent  heat  of  vaporization,  (J/kg) 

|i  =  viscosity,  (Ns/  m2) 

vj /  =  tilt  angle 
xw  =  shear  stress 


Subscripts 

c-  capillary 

e  =  electric 

h  =  horizontal  section 

h  =  hydraulic 

/  =  liquid 

max  =  maximum 

p  =  vertical  section 

v  =  vapor 
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Figure  1  The  geometry  of  the  ion-drag  pump  assisted  micro  heat  pipes 

(a)  The  schematic  of  the  system  without  electric  field 

(b)  The  look-down  configuration  of  the  system 

(c)  The  configuration  with  electric  field 
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Figure  2  (a)  Schematic  of  a  micro  heat  pipe 
(b)  Control  volume 


(a) 


Control  volume  1 
Vertical  channel 


(b) 

Control  volume  2: 
Horizontal  channel 


Figure  4  Configuration  of  the  control  volume 
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Figure  6  The  profiles  of  the  calculation  results  at  14  kV/mm,  heat  input:  41  W/cm2 
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Figure  10  Diagram  of  control  simulink  blocks 
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Figure  1 1  Temperature  response  of  control  simulink 


Figure  13  Temperature  variation  versus  time  at  15  kV/mm,  heat  flux:  40  W/crrT 
and  reference  temperature:  63  °C 
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Figure  14  Temperature  variation  versus  time  at  15  kV/mm,  heat  flux:  40  W/cm2 
and  reference  temperature:  89.3  °C. 


